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ABSTRACT 

Aims. Our goal is to determine the stellar and wind properties of seven O stars in the cluster NGC2244 and three O stars in the OB 
association MonOB2. These properties give us insight into the mass loss rates of O stars, allow us to check the validity of rotational 
mixing in massive stars, and to better understand the effects of the ionizing flux and wind mechanical energy release on the surrounding 
interstellar medium and its influence on triggered star formation. 

Methods. We collect optical and UV spectra of the target stars which are analyzed by means of atmosphere models computed with 
the code CMFGEN. The spectra of binary stars are disentangled and the components are studied separately. 

Results. All stars have an evolutionary age less than 5 million years, with the most massive stars being among the youngest. Nitrogen 
surface abundances show no clear relation with projected rotational velocities. Binaries and single stars show the same range of 
enrichment. This is attributed to the youth and/or wide separation of the binary systems in which the components have not (yet) 
experienced strong interaction. A clear trend of larger enrichment in higher luminosity objects is observed, consistent with what 
evolutionary models with rotation predict for a population of O stars at a given age. We confirm the weakness of winds in late O 
dwarfs. In general, mass loss rates derived from UV lines are lower than mass loss rates obtained from Ho-. The UV mass loss rates 
are even lower than the single line driving limit in the latest type dwarfs. These issues are discussed in the context of the structure of 
massive stars winds. The evolutionary and spectroscopic masses are in agreement above 25 Mq but the uncertainties are large. Below 
this threshold, the few late-type O stars studied here indicate that the mass discrepancy still seems to hold. 

Key words. Stars: early-type - Stars: fundamental parameters - Stars: winds, outflows - ISM: H ii regions 



1. Introduction 

Massive stars play key roles in various fields of astrophysics due 
to their strong ionizing fluxes, powerful winds and extreme stel- 
lar properties. They have been recognized as the precursors of 
long-soft gamma-ray bursts. Population III stars are thought to 
be made of very massive stars with exceptionally hard UV lu- 
minosity. As they live rapidly, they are excellent tracers of star 
formation. But they also trigger the formation of new genera- 
tions of stars through their feedback effects. The heavy elements 
produced in their cores are injected in the interstellar medium by 
stellar winds and supernovae ejecta, contributing to the chemi- 
cal evolution of galaxies. Constraining their properties and un- 
derstanding quantitatively their evolution and feedback effects is 
thus of interest for fields well beyond pure stellar physics. 

In addition to their mass, the evolution o f massive stars is 
gover ned by two m ain parameters: mass loss jChiosi & Maedeij 
and rotation dMaeder & Meynetl l2000h . The accurate de- 
termination of mass-loss rates relies mainly on two types of 
methods. First, the infrared-millimeter-radio excess emission is 
meas u red and converted to mass-loss rate ( Wright & Barlow! 
Il975t iBarlow & CohenI [19771: iLeiflierer & Robertlll99lh . This 
method suffers from little approximations but requires stars to 
be relatively nearby to be detected with the current generation 
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of instruments. Second, spectroscopic features sensitive to the 
wind extension (such as P-Cyg ni profiles or emission lines) are 
fitted by atmosphere models (IPuls et al.lll99d: iHillier & Milled 
1999). This method allows the analysis of larger samples of 
massive stars, but relies on the physics and approximations in- 
cluded in atmosphere models. The main uncertainty in mass- 
loss rate determinations (which affects both types of methods) 
is the non-h omogeneity of massive stars winds (clumping). 
Direct (E versber g et anil998.) and indirect (e.g. iHillier. ,199ll: 
iBouret et al. 2005) evidence exists for the presence of inho- 
mogeneities in O and Wolf-Ray et stars winds. Although the- 
oretically expecte d (MacGregoret al. 1979; Owocki & Rybic^ 
11984 ll985l:]Runacres & Owockil 120021) . clumping is cuiTently 
not fully und erstood and characterized, although w ork is in 
progress (e.g. iPuls et alJl2006t ISundqvist et al.ll2010l) . This is a 
severe problem for stellar evolution since mass loss is a key pa- 
rameter. Hence, having as many stars as possible studied quanti- 
tatively is desirable in order to understand their outflows. In this 
context, the late O-type main-sequence stars represent a puzzle 
of their own, since they display very weak winds with mass-loss 
rates down to a few 10"'° M© yr"' (Martins et al."2 004l2005bt 
,Marcolino et al. 2009,) . The origin of such a weakness is not un- 
derstood. 

Rotation has several effects on the appearance and evolution 
of massive stars. In addition to shaping the surface of the star 
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and modifying the surface temperature (poles being hotter than 
the equator), it triggers internal mixing processes (mainly merid- 
ional circulation and shear turbulence) which affect the transport 
of angular momentum and chemical species. Consequently, the 
timescales (such as main-sequence lifeti me) and surface abun - 
dances depend on the rotation rate (Mae der & Mevne"3l2000h . 
iHunter et alJ (l2008ll200 9) have questioned the role of rotational 
mixing in the surface chemical enrichment of B stars in the 
Magellanic Clouds and the Galaxy. They found both slowly ro- 
tating N-rich and rapidly rotating non enriched stars that, ac- 
cording to the authors, cannot be accounted for by sin gle star 
evolution with rotation. However, Maeder et al.l ([2009') argued 
that when care is taken in disentangling the various parameters 
on which the surface nitrogen content depends, 80% of the pop- 
ulation of B stars is explained by normally rotating single stars. 
Of particular importance is the age of the population: it is cru- 
cial to build samples of stars of the same age to avoid too large 
a spread in nitrogen surface abundances. So far, the relation be- 
tween chemical enrichment and rotation has not been looked at 
in more massive O stars. Since they evolve more rapidly and pro- 
duce a larger enrichment than B stars, they are important objects 
to analyze. 

The region of the Rosette nebula is especially interesting 
in this context. It harbors the ionizing cluster NGC2244 at the 
core of the nebula itself, as well as several more widely spread 
massive stars in the MonOB2 association. NGC2244 contains 
seven O stars covering the entire range of spectral types, from 
04 to 09. All of them have main-sequence luminosity classes 
(see Table [TJ. The stars of MonOB2 have similar spectral types. 
Hence, this region contains a homogeneous sample of O stars 
with different masses and luminosit ies but most lik ely of the 
same age. Binaries are also present ( Mahv et al.ll2009h . As such, 
this sample is well suited to investigate the effects of rotation 
on massive stars surface abundances. The Rosette nebula is also 
famous for the numerous young stellar objects and embedded 
cluste r s in the so-called Rosette m olecular cloud dPhelps & Ladal 
[19971 i i Francesco et al.l l2010h . Triggered star formation is 
thought to be at work in this region, due to the ionizing radiation 
and powerful winds of the NGC2244 O stars. Recent results ob- 
tained by the Herschel Space Observatory reveal the presence of 
younger objects further away fro m the HIT region center, as well 
as temperature gradients ( jd i Francesco et al.l 120101: iMotte et al] 
[20T0I ISchneider et al.ll2010ir Constraining the mechanical en- 
ergy release and ionizing flux of the massive star population is 
thus important in order to quantitatively study the process of star 
formation triggering. 

In this paper, we analyze nine O stars and one supposed B 
star in the Rosette nebula region. Seven are located in NGC2244, 
three in MonOB2. Our study encompasses six out of the seven O 
stars powering the Rosette HII regior(3. Two confirmed binaries 
(HD46149 and HD46573) and one binary candidate (HD46150) 
are included. This work comple ments and extends the analysis 
of HD48099 (M ahv et al.1 12010^, also a binary member of the 
MonOB2 association. In Sect. [2] we present the observations of 
the targets. Their analysis by means of atmosphere models is de- 
scribed in Sect. [3] The results are given in Sect. |4|and discussed 
in Sect. [5] Finally, our conclusions are drawn in Sect. [6] 



' Only the 09V star HD258691 is not included due to the lack of 
spectroscopic data. 



2. Observations and data reduction 

Table [T] summarizes the main observational properties of our 
sample stars. Photometry is take n from the Galactic O star 
Catalog of Maiz-Apellaniz et al.! (l2004h . We have adopted a 
distance of 1.55+0.15 kpc for the Mon OB2 association and 
NGC2244. Photometric dist ances range between 1 .4 and 1 .7 kpc 
according to the summary of lHensberge et al.](l2000h . We simply 
adopt the average of these values in the present work (see also 
Mahy et al. 2009). 

Spectroscopic data were obtained at the Observatoire de 
Haute-Provence (OHP, France) with the Aurelie spectrograph 
mounted on the 1.52m telescope. Several observing runs from 
2006 to 2008 allowed us to collect spectra in two different wave- 
length domains: [4450-4900] and [5500-5920]A. In the adopted 
configuration, the Aurelie spectrograph offers a resolving power 
of about = 8 000 in the blue domain whilst R = 10000 is 
reached in the red band. The exposure times, necessary to reach 
a signal-to-noise ratio higher than 150, were typically 10 to 45 
min. We also retrieved 13 spectra from the Elodie {R - 42 000) 
archives taken from November 1999 to November 2005 with the 
1.93m telescope at the OHP. The exposure times of these data 
ranged from 10 min to Ih. 

In April 2007 and March 2008, another set of spectra 
was obtained at Observatorio Astronomico Nacional of San 
Pedro Martir (SPM) in Mexico with the Espresso spectrograph 
mounted on the 2.10m telescope (R - 27 000). This echelle 
spectrograph covers the wavelength domain [3780-6950] A with 
27 orders. We used a typical exposure time from 5 to 15 min. 
Several consecutive spectra of a given night were added to im- 
prove the signal-to-noise ratio of our data. 

We have also retrieved 11 spectra from the ESO archives 
(PI: Casassus, run 076.C-0431(A) and PI: Lo Curto, run 076.C- 
0164(A)) taken with the Fiber-fed Extended Range Optical 
Spectrograph (FEROS) mounted on the ESO/MPG 2.20m tele- 
scope at La Silla (Chile). The wavelength region covered by this 
instrument is [3550-9200] A and the spectral resolution is about 
= 48 000. The data reduction was done with an impro ved ver- 
sion of the FEROS pipeline as described in lSanal (l2009h . 

Finally, to model the winds of the O-type stars with more 
accuracy, we have retrieved a total of 40 high-resolution lUE 
spectra (SWP with a dispersion of 0.2A). Such data are avail- 
able for all targets of our sample, except HD 46573. These data 
were taken between September 1978 and February 1987 with 
exposure times between 539 and 18000 s. These spectra cover 
the [1200-1800]A spectral band, which is particularly interest- 
ing for the UV P-Cygni profiles (C ivl548 - 1550 and N vl240). 

A full description of the optical instru ments used as well as 
the reduction procedure was given in Mahy et alj (l2009l) . Tables 
|2]and|3]present the main characteristics of the observ ational data. 
For th e single stars, we have selected the spectra of iMahv et al.l 
(I2009h with the widest wavelength coverage, highest spectral 
resolution and highest signal-to-noise ratio. 

3. IVIodelling 

We have used the code CMFGEN (iHiUier & MiUeii 19981) for the 
quantitative analysis of the optical and UV spectra. CMFGEN 
provides non-LTE atmosphere models including winds and 
line-blanketing. An e xhaustive description can be found in 
iHiflier & MiiieJ (Il998h and a summary of the main character- 
istics is given in iHillier et a l. (2003). CMFGEN needs as in- 
put an estimate of the hydrodynamical structure that we con- 
structed from TLUSTY models (taken from the OSTAR2002 
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Source 


ST 


U 


B 


V 


J 


H 


K 


Av 


Mv 


NGC2244 


HD 46223 


04V((f)) 


6.712 


7.484 


7.266 


6.742 


6.684 


6.693 


1.54 


-5.22 


HD 46150 


05V((f))z 


6.043 


6.875 


6.745 


6.446 


6.451 


6.453 


1.27 


-5.48 


HD 46485 


07Vn 


7.876 


8.563 


8.243 


7.565 


7.492 


7.463 


1.86 


-4.57 


HD 46056 


08Vn 


7.705 


8.445 


8.245 


7.837 


7.816 


7.837 


1.49 


-4.20 


HD 46149-1 


08V 


6.982 


7.774 


7.602 


7.245 


7.232 


7.268 


1.40 


-4.13 


HD 46149-2 


08.5-9V 














1.40 


-3.84 


HD 46202 


09.5V 


7.623 


8.360 


8.183 


7.784 


7.760 


7.737 


1.42 


-4.19 


Mon OB2 


HD 46573 


07V((f))z 


7.613 


8.273 


7.933 


7.197 


7.148 


7.145 


1.92 


-4.94 


HD 48279 


ON8.5V 


7.254 


8.047 


7.910 


7.650 


7.681 


7.710 


1.29 


-4.33 


HD 46966 


08.5IV 


5.910 


6.827 


6.877 


6.919 


6.951 


7.035 


0.71 


-4.78 


HD 48099-1 


05.5V((f)) 


5.387 


6.323 


6.366 


6.445 


6.490 


6.512 


0.73 


-5.39 


HD 48099-2 


09V 














0.73 


-3.89 



Notes. Spectral types are from lSota et al.l ( 1201 Ih and lMahv et al.l ( 120091) . Magnitudes are from lMafz-Apellaniz et al.l ( l2004f) . The extinction Ay is 
3.1 X E(B-V). The absolute magnitude assumes a distance of 1.55 kpc. For binary stars, we give the observed magnitudes of the system on the 
row of the primary star For binary systems, the absolute V magnitudes of each component are computed assuming the V-band flux ratio typical 
for stars of the same spectral type (see Sect.|4ll. 



Table 2. Journal of optical observations. 



Star 


Instrument 


Resolution 


Spectral range 


S/N 


Obs Date 


HJD 


Exposure time 








[A] 






[d] 


[s] 


HD 46056 


ESPRESSO 


18000 


[3950 - 6700] 


269 


20 Mar 2008 


2454545.6733 


2700 


HD 46149 




disentangled spectrum 


- see data in Mahv et al. (2009) 




HD 46150 


FEROS 


48000 


[4000-7100] 


310 


06 Feb 2006 


2453772.5895 


3600 


HD 46202 


FEROS 


48000 


[4000 - 7100] 


154 


05 Jan 2006 


2453740.6191 


600 


HD 46223 


FEROS 


48000 


[4000-7100] 


194 


08 Feb 2006 


2453774.5460 


3600 


HD 46485 


FEROS 


48000 


[4000-7100] 


166 


05 Jan 2006 


2453740.6290 


600 


HD 46573 


FEROS 


48000 


[4000-7100] 


268 


07 Feb 2006 


2453773.5904 


3600 


HD 46966 


FEROS 


48000 


[4000-7100] 


200 


03 Jan 2006 


2453738.7201 


600 


HD 48099 






disentangled spectrum - 


see Mahv et al. 


(2010) 




HD 48279 


ELODIE 


42000 


[4000 - 6800] 


214 


25 Jan 2000 


2451569.5028 


2700 



Table 3. List of lUE UV spectra. 



Star 


data ID 


Obs. Date 


HJD 


Exposure time 








[d] 


[s] 


HD 46056 


SWP06949 


22 Oct 1979 


2444168.92297 


5100 


HD46149 


SWP28151 


11 Apr 1986 


2446532.40344 


2640 


HD46150 


SWP10758 


05 Dec 1980 


2444578.86527 


1200 


HD 46202 


SWP30299 


12 Feb 1987 


2446839.37549 


5100 


HD 46223 


SWP10757 


05 Dec 1980 


2444578.81523 


1980 


HD 46485 


SWP28196 


19 Apr 1986 


2446539.98248 


10800 


HD 46966 


SWP24196 


18 Oct 1984 


2445991.75211 


540 




SWP27918 


15 Mar 1986 


2446505.35266 


540 


HD 48099 


disentangled spectrum - see Mahv et al. (2010) 


HD 48279 


SWP06504 


14 Sep 1979 


2444130.71788 


7200 



Notes. lUE UV spectra retrieved from the archive and used for the analysis. The wavelength range is 1 150-1975 A and the spectral resolution 0.2 
A (/f=7500 at 1500 A) for all spectra. For HD 48279, we used the average of the two spectra listed in the table since they were obtained in the 
same mode (high aperture), with the same exposure time and do not show any variability. For the other stars, we chose the spectra with the highest 
exposure time in the high aperture mode. 



grid of iLanz & Hubenvll2003h connected to a y6 velocity law of 
the form v - Voo(l - R/r'f where Voo is the wind terminal ve- 
locity. We ado pted /? - 0.8 since t his is the typical value for 
O dwarfs (e.g. iRepolust et al]|2004l) . Once the model was rea- 
sonably converged, the radiative acceleration was computed and 
used to iterate on the initial hydrodynamical structure in the in- 
ner atmosphere. The radiation field and level populations were 



then converged once again. Two to five of these global hydrody- 
namical iterations were performed before the atmosphere model 
finally converged. 

Our final models include H, He, C, N, O, Ne, Mg, Si, S, 
Fe with the solar composition of iGrevesse et al . (2007) unless 
otherwise stated. The super-level approach is used to reduce the 
amount of memory requirements. On average, we include 1600 
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super levels for a total of 8000 levels. For the formal solution of 
the radiative transfer equation leading to the emergent spectrum, 
a microturbulent velocity varying linearly (with velocity) from 
10 km s"' to 0. 1 X Voo was used. 

We included X-ray emission in the wind since this can af- 
fect the ionization balance and the strength of key UV diagnos- 
tic lines. In practice, we adopted a temperature of 3 million de- 
grees and adjusted the flux level so that the X-ray flux eff'ectively 
coming out of the atmosphere matches the observed Lx/L\,o\ ra- 
tio. This ratio was computed from the observed X-ray fluxes of 
IWang et al. (2008) and our derived luminosities. In case no X- 
ray flux was d etected, we simply adopted the canonical value 
Lx/Lboi = 10"^ (ISana et al.l |2006; Naze 2009) . 

In practice we proceeded as follows to derive the stellar and 
wind parameters: 

• Effective temperature: we used the classical ratio of Hei to 
Hen lines to constrain T^g. The main indicators were He: 
4471 and Hen 4542. Additional diagnostics are Hei 4026 
(blended with Hen 4026), Hei 4388, Hei 4713, Hei 4921, 
He I 5876, He it 4200, He n 5412. The typical uncertainty on 
the Teft determination is 1000 K. 

• Gravity: The wings of the Balmer lines Hy6, Hy and Wd are 
the main \ogg indicators. An accuracy of about 0.1 dex on 
\ogg is achieved. 

• Luminosity: We fitted the flux calibrated lUE spectra and 
UBVJHK fluxes to derive both the luminosity and the ex- 
tinction. For this, we adopted a distance of 1.55+0.15 kpc 
(see dis cussion inlMa hy et al. 2009h . The Galactic extinction 
law of iSeatonl (Il979h and iHowarthl(ll983h was used. We ad- 
justed the stellar luminosity and E(B-V) to reproduce the 
UV-optical-infrared SED. In those cases where no UV spec- 
tra were available, we simply derived the extinction from 
E(B-V) and adjusted the luminosity to reproduce the V mag- 
nitude. 

• Wind terminal velocity: the blueward extension of P-Cygni 
lines observed in the UV provides Voo with an accuracy of 
100 km s '. 

• Mass-loss rate: We used two diagnostic lines to constrain the 
mass loss rate: the UV P-Cygni profiles and Ha. In principle, 
a single value of M should allow a good fit of both types of 
lines. However in practice, several recent studies have shown 
that this was not the case. In order to be as complete as pos- 
sible, we thus provide both values: the "UV mass loss rate" 
(Muv, given in TableHJ and the "Ha mass loss rates" (Maa, 
given in Table |5]l. 

The projected rotational velocities were derived by direct 
comparison of our synthetic spectra to observed line profiles for 
fast rotating stars (i.e., those with V sin/ > 150 km s '). The 
rotationally broadened synthetic line profiles usually provided a 
good match to the observed profile. In the case of moderately 
to slowly rotating objects, some amount of macroturbulence had 
to be introduced to coiTectly reproduce the line profiles of sev- 
eral features (e.g., Hei 4713, Civ 5812, Hei 5876). The need 
for extra broadening is well documented (Howarth et al. '19971 
Ryans et al. 2002; Howarth et al. 2007; Nieva & Przybilla 20ol 
Martins et al.ll2010t IFraser et al.ll2010 ) but its origin is unclear. 
Aerts et al.l (l2009h recently suggested that non-radial pulsations 
could trigger large-scale motions (hence macroturbulence), but 
this needs to be confirmed by a study covering a wider param- 
eter space. A recent study by Simon-Diaz et al. (2010) showed 
that the amplitude of macroturbulence was correlated to the am- 
plitude of line profile variability in a sample of OB supergiants. 



In practice, we introduced macroturbulence by convolving our 
(rotationally broadened) synthetic spectra with a Gaussian pro- 

file (oc e -'moc where v^ac is the macroturbulent velocity), thus 
mimicking isotropic turbulence. This is obviously a very sim- 
ple approach, but it significantly improves the quality of the fits. 
Given the present limitations, we restricted ourselves to rough 
estimates of the amount of macroturbulence by judging the fit 
quality by eye. In practice, we used He 1 47 1 3 as the main indica- 
tor of macroturbulence since it is present with suflicient SNR in 
all our sample stars. Secondary indicators were the C iv doublet 
at 5801-5812 A and He: 587 6. V sin; was also obtained from 
the Fourier transform method dSimon-Diaz & H errero 2007). In 
practice, for the moderately to slowly rotating stars, we obtained 
upper limits from the absence of zero and the position of the 
noise level in the FT. 

We also determined the nitrogen content of our sample stars. 
We relied mainly on the N iii lines between 4500 and 4520 A. 
They are present in absorption in all the stars. They are not af- 
fected by winds and are strong enough for abundance determi- 
nation. The uncertainties are of the order of 50%. They were 
estimated from the comparison of the selected N iii lines to mod- 
els with various N content. The errors do not take into account 
any systematics related to atomic data. 

When possible, the degree of inhomogeneities of the winds 
of our sample stars was determined. Clumping is implemented 
in CMFGEN by means of a volume filling factor / following the 
law / = /oo + (1 - fcc)e^^'^^'''' where fco is the maximum clump- 
ing factor at the top of the atmosphere and Vd a parameter in- 
dicating the position w here the w ind starts to be significantly 
clumped. As shown by i Bouret et al. (2005), Ov 1371 and Niv 
1720 are two UV features especially sensitive to wind inhomo- 
geneities in early O stars. We have used these lines to constrain 
fco in the earliest O stars of our sample. In mid to late O stars, no 
UV clumping diagnostic has been identified. This is confirmed 
by our study. For the 07-09 stars of our sample, we computed 
models with several values of /«, (scaling M so that M/ ^ff is 
constant) and found no difference in the UV wind sensitive lines. 
Hence, we decided to adopt /oo = 1 (i.e. homogeneous model). 
Our mass-loss rates for those objects should thus be regarded as 
upper limits. 

4. Results 

The derived stellar and wind parameters of the sample stars are 
gathered in Table |4] The values of the Ha mass loss rates are 
given in Table |5] We show the best optical fits for HD 46223 
and HD 48279 in Figs. [1] The optical fits for the other stars are 
gathered in Appendix A. The fits to the UV spectra for all stars 
are shown in Fig. [2] Below, we briefly comment on each star 

HD 46223: a clumped wind (/=0.1) is required to fit Ov 
1371. The nitrogen content is difficult to estimate given the 
weakness of the N iii 4500^520 features. A value of N/H-1.0 
X 10""* is preferred. It also provides a good fit to the Niv 1720 
feature. The parameters are similar to those of iMartins et al.l 
(l2005bl) within the uncertainties. The better quality and wider 
wavelength coverage of the present optical spectra explains the 
refined effective temperature estimate. 

HD 46150: this is a binary candidate. We thus derive upper 
limits on the luminosity and mass loss rate. A clumping factor 
of 0.1 is necessary to correctly reproduce the Ov 1371 A line 
profile. Models with higher (lower) clumping factor yield a too 
strong (weak) absorption in the blue part of the profile. See also 
Sect. l5.2l for further discussion. 
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Table 4. Derived stellar properties of the sample stars. 



Source 


ST 


7'eff 




log(g) 


log(Muv) 


u 




V sin; 








N/H 


log(Qo) 






[kK] 








[kms-'] 


[km s-'] 


[km s-'] 


[Mq] 


[Mq] 


[xlO-t] 




NGC2244 


HD 46223 


04V((f)) 


43.0 


5.60+0.11 


4.01 


-7.17 


0.1 


2800 


100 


32 


52 1+^- 


48.3 + 19.3 


7.0+2.0 


49.40 


HD 46150 


05V((f))z 


42.0 


<5.65 


4.01 


-7.30 


0.1 


2800 


100 


37 


< 53.6 


59.5 + 23.7 


3.0+2.0 


49.36 


HD 46485 


07Vn 


36.0 


5.05+0.11 


3.85 


-7.80 


1.0 


1850 


300 




28.6+" 


19.5 + 7.9 


<1.2 


48.62 


HD 46056 


08Vn 


34.5 


4.85±0.12 


3.89 


-8.50 


1.0 


1500 


330 




23. 1«-^ 


15.8 + 6.7 


<0.6 


48.32 


HD 46149-1 


08V 


37.0 


4.90±0.12 


4.25 


< -9.0 




1300 


~0 


24 


25.9+3-^ 
20.6:|* 


30.9 ±21.6 


0.8+0.5 


48.42 


HD 46149-2 


08.5-9V 


35.0 


4.65+0.16 


4.01 








100 


27 


12.5 + 10.1 


0.8+0.5 


48.10 


HD 46202 


09.5V 


33.5 


4.85+0.12 


4.10 


-9.0 


1.0 


1200 


20 


17 


99 7+2.0 
^^•'-2.4 


29.0 + 12.4 


1.0+0.5 


48.19 


Mon 0B2 


HD 46573 


07V((f))z 


36.5 


<5.20 


3.75 








20 


43 


< 32.5 


20.6 + 8.3 


3.0+1.0 


48.84 


HD 48279 


ON8.5V 


34.5 


4.95+0.11 


3.77 


-8.8 


1.0 


1300 


125 


22 


24.4+2-5 


15.2 + 6.1 


5.0+3.0 


48.47 


HD 46966 


08.5IV 


35.0 


5.20+0.11 


3.75 


-8.0 


1.0 


2300 


50 


27 


55 3+^-^^ 

-iO 


24.5 + 9.9 


1.2+0.5 


48.70 


HD 48099-1° 


05.5V((f)) 


44.0 


5.65+0.07 


4.50 


-7.60 


0.01 


2800 


91 





155.0 + 98.7 


5.0+2.5 


49.36 


HD 48099-2° 


09V 


32.0 


4.60+0.06 


3.51 








51 







5.0 + 3.5 


0.6+0.3 


47.93 



Notes. The columns give: name, spectral type, effective temperature (uncertainty lOOOK except for HD 46149-1 / HD 48099-1, 2000K and HD 46149-2 / HD 48099-2, 3000K), luminosity, 
effective gravity (uncertainty 0.1 dex, except for HD 46149 and HD 48099, 0.2 dex), UV mass-loss rate, terminal velocity, projected rotational velocity, macroturbulent velocity, evolutionary mass, 
spectro scopic mass, nitrog en content, ionizing flux. Spectroscopic masses are computed from the true gravity (i.e. log g corrected from the effect of centrifugal force). 
from lMahv et alj ( 1201 Oh . except that the error bars on Teff and logg are similar to HD 46149 (see text for discussion). 
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Fig. 1. Best CMFGEN fits (red) of the optical spectra (black) of stars HD46223 (top), HD48279 (bottom). The models have the 
mass loss rate best reproducing the UV features. The unfitted feature around 4425 A is a DIB. 
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Fig. 2. Best CMFGEN fits (red) of the UV spectra (black) of the sample stars. 



HD 46485: this is the second fastest rotator of our sample, 
with V sin/ - 300 km s The fit of the optical and UV spec- 
tra is of excellent quality. The only notable exceptions are the 
cm line complex at 4640 A, Civ 5801-5812 A and He i 5876 
A. The C m 4650 line complex is slightly too weak in our best 
fit model. However, we found out that this transition is sensi- 
tive to detailed line-blanketing eff'ects. The lower level of these 
transitions is connected to the C iii ground level by a transition 
at 538A. This transition is blended with Fe iv transitions. The 
interaction between these UV lines thus aff'ects the strength of 



the Cm 4650 lines , in a manner similar to that described by 
iNaiarro et al.l (l2006h for the He i singlet lines. 

HD 46056: this is the fastest rotator of our sample and no 
macro turbulence is necessary to reproduce the line profiles. The 
fit is good with the exception of the Cm 4650 line complex 
which is underpredicted. The Cm 4650 lines suffer from the 
same problem as for HD 46485. 

HD 46149: this is a binary (iMahv et al.ll2009h . The spec- 
tra have b e en di sentangled using the method described in 
Ma hv et al.l (l2010l) . They have subsequently been used to clas- 
sify the two components of the system. The classical diag- 
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nostic log ^w(Hen4542) ^^^'^ ^ spectral type 08 for the pri- 
mary and 08.5-09 for the secondary. Similarly, a luminosity 
class V was assigned based on the values of log 
(IConti&Alschuleiill971h . These classifications are slightly dif- 
ferent from those derived by Mahy et al (2009) although for 
the primary it remains in the error range. However, such a dis- 
crepancy for the secondary could come from the disentangling 
process, making the disentangled spectra uncertain. Indeed, they 
were computed from a sample of observed spectra taken essen- 
tially at similar phases where the separation between the radial 
velocities of both components (RVs - RVp) are of about -30, 
and -h140 km s Moreover, the orientation and the large ec- 
centricity of the orbit imply an asymmetric excursion in radial 
velocity which could also affect the accuracy of the broadest 
lines. We thus caution that the disentangling process does not 
provide perfect results and should be refined in the future when 
more data taken at different phases will be obtained. This is es- 
pecially true for the wings of broad lines which need to be cor- 
rectly sampled and fully separated in the observed spectrum (at 
maximum separation) to be reconstructed coiTectly (see below - 
HD 48099). Consequently, the gravity and spectroscopic masses 
in Table |4] are highly uncertain. The effective temperature and 
nitrogen content have also larger error bars than the other stars 
of our sample: 2000 K (3000 K) on Teff for the primary (sec- 
ondary) and ~70% on N/H. In addition, given the poor quality of 
the spectrum o f the secondary, w e adopted log g = 4.0. Using the 
calibrations of iMartins & Pled (j2006), we find that the V-band 
flux ratio between a 08V and a 09V star is x 1.3. Optical pho- 
tometry of HD 46149 indicates E(B-V)= 0.45, leading to Ay = 
1 .40. For the adopted distance, the absolute V-band magnitude 
of the system is thus My = -4.75. Given the V-band flux ratio 
of the two components, we have My™"'^ - -4.13 and ^^'='=°"<''"t 
= -3.84. 

HD 46202: the parameters derived f or HD 46202 are very 
similar to those of iMartins et alj (l2005bh . The fit to the optical 
spectrum is of excellent quality. A combination of a low rotation 
rate (20 km s"' ) and a macroturbulent velocity of 17 km s"' cor- 
rectly reproduces the photospheric line profiles. See also Sect. 
l5.2l for further discussion. 

HD 46573: this is a binary candidate. No UV spectra are 
available for that star, so that our analysis is restricted to the op- 
tical range. A very good fit is achieved. Only the core of the 
Balmer lines is too deep in our synthetic spectra. The line pro- 
files show that macroturbulence is important. We found that a 
combination of V sin/ in the range 40-60 km s"' with a macro- 
turbulent velocity of 40 to 50 km s"' leads to much better fits 
compared to purely rotationally broadened profiles. 

HD 48279: our synthetic spectrum provides a good fit of 
most optical and UV lines. The singlet lines Hei 4388 A and 
He 1 492 1 A are too weak in our spectrum. They are know n to be 
sensitive to details of line blanketing (Nai arro et alj|2006h . 

HD 46966: the fit of the optical and UV spectra is one of the 
best of the present work. The only problem is observed in He n 
4686 A for which our synthetic profile is too deep. A projected 
rotational velocity of the order 40-50 km s"^ combined with a 
macroturbulent velocity ~ 30^0 km s"' provided better fits than 
models with only rotationally broadened profiles. 

HD 480 99: this star is a binary and has been studied by 
iMahv et alJ (golO). We have not re-analyzed the two compo- 
nents and we have simply adopted the parameters derived by 
Mahy et al. However, we add a few words of caution. The spectra 
used in the disentangling process of this binary system sample 
the full orbital cycle and hence do not suffer from the same lim- 



itations as in the case of HD46149. But whilst we expect narrow 
spectral features to be reliably reconstructed in the disentangled 
spectra of both the primary and secondary star, the reconstruc- 
tion of the broad features is more uncertain . Indeed, given th e 
rather low orbital inclination of HD48099 dMahv et al.ll2010l) . 
the maximum radial velocity excursions of the stars remain small 
(RV'^"' - RV'^'"' ~ 155 km s"' - RV'^"' and RV^"" being the 
maximum radial velocities of the primary and secondary com- 
ponents) compared to the possible widths of broad Balmer lines 
(±950km s"' for Hy). This poor sampling of the broadest lines 
can lead to a poor reconstruction of the wings of the Balmer 
lines especially for the fainter secondary component. Hence, the 
gravity of the latter is very likely underestimated. There is un- 
fortunately no way to remedy to this situation and as a result, the 
spectroscopic masses for this system are unreliable as becomes 
clear for instance when comparing the spectroscopic mass ra- 
tio from Table 4 (M1/M2 ~ 30) with the dynamical mass ratio 
(Ml /M2 - 1 .77) from the orbital solution of Mahy et al.. (2010.) . 
In view of this discussion, we chose to adopt the same error bars 
on Teff and logg as for HD 46149 (i.e. 2000 K on T^ff - 3000 K 
for the second ary- and 0.2 dexon log g). The error bars are thus 
larger than in iMahv et al.l (l2010l) . They are also realistic, since 
reducing logg by 0.25 dex, we find that we need an effective 
temperature lower by 1000-2000 K to correctly reproduce the 
He I and He 11 lines of the primary. 



5. Discussion 

5.1. HR diagram and evolutionary status 

The HR diagram of the stars analyzed in the present study is 
shown in Fig. [3] Red/pink (blue/cyan) colors are used for the 
stars of NGC2244 (Mon OB2 respectively). 

The stars of NGC 2244 have an age less than 5 Myr 
The most massive stars of the cluster, namely HD 46223 and 
HD 46150, are less than 2 Myr old. There is no sig nificant age 
difference between these two stars, contrary to what Wang et al.l 
(2008) suspected. Whether the trend for the most massive stars 
being younger than the lower mass O stars is real or an artifact 
of the parameter determination is unclear. It is known that the 
stellar parameters of an average early O star indic ate a younger 
age than that of late-type O stars (e.g. Fig. 12 of IMartins et al.l 
2005a). An explanation to this trend is that late-type O stars have 
lower ionizing fluxes and weaker winds than early-type O stars. 
Consequently, they remain hidden in their parental cloud for a 
longer time, so that on average late-type O stars appear older In 
the present case, since both early and late-type stars are visible 
at the same time, one needs another explanation. A possibility is 
that the most massive stars were the last to form. Alternatively, 
there might be a bias in the determination of the effective tem- 
perature of the hottest O stars. It would have to be 2000-3000 K 
lower in order to lead to an age similar to the late-type O stars. 
Such a reduction is difficult to explain. Line-blanketing, when 
included in model atmospheres, leads to a downward revision of 
effective temperature s by several thousands K. But as shown by 
Martins etall (l2002h . C, N, O, Si, S and Fe are the main con- 
tributors to opacities in O stars atmospheres. Since our present 
models include all these elements and additional ones, it is un- 
likely that we overestimate Tetj by 2000-3000 K. Since we are 
using clumped models for the hottest stars and homogeneous 
winds for the coolest ones, one might wonder whether this af- 
fects the temperature determination. However the photospheric 
lines used to constrain Teff are unaffected by clumping since they 
form deep in the atmosphere where the wind is still homoge- 
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Fig. 3. HR diagram of NGC2244. The ev olutionary tracks 
and isochrones are from iMevnet & MaedeJ |2003). They are 
computed for an initial rotational velocity of 300 km s"^ 
Triangles are the components of the binaries HD48099 (cyan) 
and HD46149 (pink). Stars belonging to NGC2244 are shown in 
red/pink and stars from the Mon OB2 association are shown in 
blue/cyan. Aitows indicate upper limits on luminosities for SBl 
and binary candidates. 



neous even in clumped models. Another speculative hypothesis 
is that the most massive O stars have, on average, larger rotation 
rates. If so, the use of the same set of evolutionary tracks (i.e., 
computed with the same initial rotational velocity) for early and 
late-type O stars is inappropriate. Indeed for faster rotation, the 
main-sequence lifetime is increased. iMaeder & Mevned (I2OO0I) 
concluded that the use of non-rotating tracks could le ad to an un- 
deres timate of the actual age of stars by about 25%. lBrott et al.l 
(1201 lai) show the effect of rotation on isochrones (see their Fig. 
7) and confirm fast rotating massive stars stay closer to the zero- 
age main sequence compared to moderately rotating lower mass 
stars of the same age. This could partly explain the age differ- 
ence we observe between early and late O stars. But this relies 
on the assumption that rotation is faster for larger masses among 
O stars. So far, there is no detailed study of the evolution of ro- 
tational velocities with spectral type among O sta rs. An averag e 
V sin/ of 129+82 km s ' results from the work oflPennyl(ll996l) . 
consistent with the more recent result of iPennv et al.l (12004 ) - 



131+93 km s-\ iHowarth et all (11997 ) found a distrib ution of 
V sin/ peaked at about 100 km s~' . Concerning B stars. . Abt et al.l 
(l2002h found values of 100-130 km s^'for B0-B5V stars. It is 
thus not clear whether there is a trend of larger rotation rate be- 
tween B and O stars (and consequently not clear whether such a 
trend continues within the O class). 

The stars in the Mon OB2 association appear to have an 
age of about 1-5 Myr, with a prefeiTed value between 2 and 
4 Myr The only exception is the primary of the binary system 
HD 48099. As akeady discussed by Mahy et al. (2010), several 
factors can explain the apparent younger age of this star: rapid 
rotation, mass and angular momentum transfer, inappropriate 



evolutionary tracks. Rapid rotation was the favored explanation. 
The present study confirms that the secondary star has an age 
similar to the other stars of the association, making the primary 
peculiar It is thus very likely that the position of the primary in 
the HR diagram is due to a rapid rotation and that comparison 
with normal (i.e. moderate rotation) evolutionary tracks is not 
appropriate. 

5.2. Wind properties 

In Table|5] we gather the mass loss rates derived from our quan- 
titative analy sis (M\w a n d Mun ) as well as the theoretical mass 
loss rates of IVink et al.l (1200 lb and the values L/c^. The latter 
corresponds to the mass loss rates obtained if the driving is due 
to a single line located at the emission peak of the SED. They are 
thus a lower limit to the mass loss rate expected by the radiation- 
driven wind theory. From Table |5] we can draw the following 
conclusions: 

1. All stars have Muv < Mho- The difference is between a fac- 
tor of 10 and a factor of a few 100! 

2. The theoretical mass loss rates of IVink et al.l (1200 Ih are sys- 
tematically larger than the UV mass loss rates. The discrep- 
ancy reaches 2 orders of magiiitudes in the latest type O stars , 
as usually seen (iMartins et al.! 2005bVM arcolino et al.l2009l) . 
The Ha mass loss rates are usually in rather good agreement 
with the theoretical predictions (when comparisons are made 
for unclumped M for the early O dwarfs). 

3. For all late type O dwarfs, the UV mass loss rates are lower 
than the single line driving limit. This is a serious puzzle. 
If the UV mass loss rates are the correct ones, this implies 
that line driving is less efficient than thought. It could mean 
that the line absorbs only at discrete values of velocities in 
the range Q-Vca (i.e., that absorbing material is present only 
at specific velocities, see below). For the earliest O dwarfs, 
the UV mass loss rates are only a factor 3-4 larger than this 
limit. 

The results of Table |5] raise the question of the validity of 
our mass loss rate determinations. We list a few reasons why 
they could be biased: 

- The mass-loss rates are based on a few features, and some- 
times based on only a single line (although the absence of 
other wind features provides upper limits). 

- The inferred mass-loss requires that the wind ionization 
structure derived by CMFGEN is accurate. Since CMFGEN 
provides a good fit to the photospheric spectrum, there is 
no a priori reason that CMFGEN should get the ionization 
structure dramatically wrong. Studies of Hii regions gener- 
ally find that the far UV spectrum predicted by CMFGEN is 
globally consistent with that required to ionize the n ebul ae 
(iMorisset et ani2004t ISimon-Dfaz & Sta siiiska 2008). 

- CMFGEN does not compute the wind properties from first 
principles (see Sect. |3]l. The mass loss rate and velocity 
law are adopted as input to compute the atmosphere struc- 
ture and the emerging spectra. However, the results can be 
used to check the consistency between the adopted mass- 
loss rate/velocity law and the momentum absorbed by the 
wind (the knowledge of the atmospheric structure allows the 
computation of the radiative forces). 

- Recent developments in the understanding of the clumping 
properties of O stars have shown that porosity (and its ef- 
fect on the velocity structure, called vorosity) could alter 
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Table 5. Summary of various mass loss rates for our sample stars. 



Star 


ST 


log(Muv) 


log(4) 


log(MH„) 


\og(Mvi„k) 


HD 46223 


04V((f)) 


-7.17 


-7.57 


-6.20 


-5.75 


HD 46150 


05V((f))z 


-7.30 


-7.52 


-6.40 


-5.68 


HD 46485 


07Vn 


-7.80 


-8.12 


-6.45 


-6.50 


HD 46056 


08Vn 


-8.50 


-8.32 




-6.80 


HD 46149-1 


08V 


< -9.0 


-8.27 




-6.74 


HD 46202 


09.5V 


-9.00 


-8.32 


-7.10 


-6.72 


HD 46573 


07V((f))z 




-7.97 


-6.30 


-6.25 


HD 48279 


ON8.5V 


-8.80 


-8.22 


-6.80 


-6.64 


HD 46966 


08.51V 


-8.00 


-7.97 


-6.40 


-6.45 


HD 48099-1 


05.5V((f)) 


-7.60 


-7.52 




-5.66 



Notes. The columns give: name; spectral type; derived UV mass loss rate; mass loss rate expect ed for driving by a single line located at the SED 
emission peak in the single scattering condition; Ho- mass loss rate; theoretical mass loss rates o f IVink etalJ ( l200lh . The filling factor /„ (see Sect. 
O is 0. 1 for HD 46223 and HD 46150, 0.01 for HD 48099-1, and 1.0 (i.e. no clumping included) for all other stars. 



the shape of key diagnostic Unes (e.g. lSundq vist et a02010h . 
Neglecting this effect could thus bias our determination. 

In the following we investigate in more detail the last two 
points. To check the hydrodynamics of our models, we have 
run several tests on HD 46202 which is typical of the late type 
O stars with only weak evidence for winds in Civ 1548-1550 
(and perhaps N v 1240). We derive a UV mass-loss rate of 1.0 
xIO^^Mq yr"'. This is a factor of 5 lower than the single-line 
limit (4.8 xlO^'^M© yr"'), and a factor of 190 lower than that 
derived using Vink's prescription for estimating the mass-loss 
rate (which gives M=1.9 xlO"^ Mq yr"'), and is primarily based 
on the strength of the absorption associated with the Civ 1548- 
1550 doublet. As noted previously, the single line limit is larger 
than the derived UV mass loss rate. To check the hydrodynam- 
ics, we ran additional models with as many elements/lines as 
allowed by our computational resources (the models include H, 
He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni and a total of about 
10 000 atomic levels) to be able to calculate the radiative force 
as accurately as possible. 

The momentum equation can be written as: 

dV -1 dP 

v-r^—^-s + Sr (1) 

dr p dr 
where 

g = GMIr" (2) 
and, gr, the radiative acceleration, is given by 
4;r r 

gr = —P I XvHydv. (3) 

Once gr is computed, we can thus check whether Eq.[T]is satis- 
fied or not. 

When clumping is important we utilize the same equation, 
however this is only an approximation. First we note that the 
presence of clumps implies a dynamic flow, but the above equa- 
tion is for a steady flow. Second, it is unclear how to handle the 
pressure term when clumping is important - in practice this is 
not crucial for O stars since we assume the clumping to be ini- 
tiated close to, or above the sonic point where the gas pressure 
term can be neglected. The expression for gr, p and;t'y are eval- 
uated in the clumps. 



In our best fit model to the UV lines, the radiation field ex- 
ceeds that required to drive the wind everywhere - at most lo- 
cations by over an order of magnitude. This result is not unex- 
pected, since our derived mass-loss rate is lower than the single- 
line limit. As a check of the consistency of the calculations we 
computed a model using a mass loss rate of 4.0 xlO"^ Mq yr"'- 
this is a factor five below Vink's estimate, but a factor of 40 
above what we derived. As shown in Fig. |4] (right panel), the 
predicted UV spectrum shows some crucial disagreements with 
observation. Particularly noteworthy is the presence of a satu- 
rated Civ profile, and two high-velocity wind features due to 
Si IV. In the observed spectrum. Si iv shows only photospheric 
component^ The high mass loss rate model shown has no X- 
rays — when we include X-rays Nv 1240, unlike the observa- 
tions, shows a strong P Cygni profile. The weakness of the N v 
profile relative to C iv is consistent with the late spectral class of 
HD 46202, and argues against a shift of the wind to high ioniza- 
tion stages. 

For a mass-loss rate of 4 xlO"** Mq yr"', the wind momen- 
tum deposited is close to that needed to drive the wind for veloc- 
ities in excess of 700 km s"', and for a standard P - I velocity 
law. Below 700 km s"' there is insufficient line force, sometimes 
by a large factor, to drive the flow. This is an important issue 

— for a steady flow the momentum balance must be satisfied 
at all depths, not just globally. Clearly we have an inconsistency. 
Either the derived UV mass-loss rate is wrong or there is an addi- 
tional retarding force that needs to be included in the momentum 
balance equation. 

For early type O dwarfs, mass loss rates are larger and closer 
to the theoretical predictions. However, problems still occur. We 
take HD 46150 as an illustrative case. For that star, we derive 
a mass-loss of 5.0 x 10 Mq yr ' from UV diagnostics. This 
gives a Civ profile in reasonable agreement with observation, 
but the N v doublet is overestimated. Surprisingly both the C iv 
and N v profile do not reach zero intensity (C iv has a residual 
intensity of 15% while the residual intensity for N v is 35%). 

On the other hand, the mass-loss derived from Ha (assum- 
ing /=0. 1, Vc/=100 km/s) is nearly an order of magnitude larger 

- 4.0 x 10"^ Mq yr"'. This mass-loss rate also provides a bet- 
ter fit to the cores of H/3 and H6. The optically derived mass- 
loss is inconsistent with the O iv 1339, 1 343 and O v 1371 lines - 
the theoretical spectra show noticeable perturbations by the wind 
whereas the observed profiles are photospheric only (see Fig.|4]l. 



^ If we adopted a lower T^if for HD 46202, the disagreement with 
Si IV would be even more striking 
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The fit to C iv is somewhat worse although definitive statements 
are impossible since C iv shows a strong P Cygni profile that is 
insensitive to M. 

Clumping, and more precisely porosity, might be a way to 
partly solve the issues mentioned above. Porosity comes in two 
flavors — spatial and velocity. Spatial porosity has been invoked 
to help explain the observed X-ray profiles, which tend to show 
much more symmetric pro files than expected on the basis of 
Ha derived mass-loss rates (lOskinova et al.ll2b07h . An alterna- 
tive explanatio n is that Ha m ass-loss rates are too high — by a 
factor of a few (ICohen et al. t2010). T he need for spatial po rosity, 
and its importance, is controversial. ICohen et al.l (|2010) argue 
that the the X-ray profiles in ^ Pup are much better explained by 
a simp le reduction in mass-loss. In addition, lOwocki & CohenI 
(I2006h have argued that the wind clumping structure is inconsis- 
tent with that required to give rise to a significant porosity effect. 

Vorosity (i.e. porosity in velocity space) can arise in a 
clumped medium because the clump velocities are discrete - 
along a given sight line not all velocities may be present. Since 
continuum transfer is unaffected by the presence of velocity 
fields, vorosity, unl ike spatial poro s ity, on ly affects line pro- 
files|3 As shown bv iPrinia & Massal (l2010l) . vorosity could ex- 
plain that the ratio of optical depth in the two components of 
Si IV 1393,1403 in B supergiants is not equal to the ratio of os- 
cillator strengths, as expect ed in the case of h omogeneous flows. 

Recent simulations by ISundqvist et al.l (12010. .201-1.) have 
highlighted the importance of vorosity (among others) on the 
formation of strong UV lines. They show that accounting for 
vorosity, UV line profiles can be strongly reduced for a given 
mass loss rate. At the same time. Ha is less affected. This tends 
to bring the optical and UV mass loss rates into better agree- 
ment when fitting oberved spectra. But these results also indi- 
cate that if vorosity is important, it implies that the winds are 
clumped. Consequently the Ha mass loss rates derived using 
homogeneous models (such as most values in Table |5]l are up- 
per limits. These computations are still exploratory (there is no 
detailed radiative transfer nor feedback on the atmosphere struc- 
ture) but they appear to be a way to solve the issues presented 
above. 

Aside from this, the effect of porosity and vorosity on the- 
oretical mass-loss rates is unclear. If it begins beyond the crit- 
ical point we would expect a reduction in the line force, and 
hence a reduction in the terminal velocity of the wind, with no 
change in the mass-loss rate. The reduction in line force occurs 
because thick lines would absorb less continuum flux, due to 
the presence of absorbing material only at specific velocities and 
not all over the entire 0-Voo velocity range. The radiation force 
from thin lines would not alter since only the column density, 
not the velocity field, is important. Of course the definition of 
which Unes are thin and which are optically thick may be af- 
fected by the wind structure. If vorosity is important below the 
critical point we might expect a reduction in the mass-loss rate. 

Clearly, the structure, properties and dynamics of O star 
winds is complex and requires improvements in both spectro- 
scopic analysis and theoretical predictions. But vorosity appears 
as a valid way of solving (at least partly) some of the issues 
raised above. 

Whatever the exact nature of these winds, some qualita- 
tive features can be obtained directly from the observations. 
This is what we show in Fig. |5] Most of the results presented 
above rely directly or indirectly on atmosphere models. In Fig. 



' For the X-ray lines vorosity is not important - it is the continuum 
opacity which primarily affects the shape of these lines. 



|5] we show a direct comparison between the observed spectra 
of two stars: HD46966 and HD48279. According to Table H 
they have very similar Tetj and gravities. This is mirrored in 
Fig. |5]by the perfect match of the He: and Hen features and 
the wings of Balmer lines. However, Nv 1240 and Civ 1548- 
1550 show a smaller blueward extent in HD48279 indicating a 
lower terminal velocity. In addition. Civ 1548-1550 is clearly 
weaker in HD48279 while the Nv 1240 profile are of similar 
strength (slightly weaker in HD48279 though). The Ha pro- 
files are barely differenciated. From the C iii optical features, we 
tentatively estimate the carbon abundance (C/H) of HD48279 
(respectively HD 46966) to be 1.0 x lO"'* (resp. 2.0 x lO""*). 
Consequently, HD 48279 exhibits a larger N/C ratio compared 
to HD46966. If the mass loss rate were the same in both stars, 
one would expect a weaker C iv 1548-1550 line but also a much 
stronger N v 1240 line (the profile is not saturated) since the ni- 
trogen overabundance in HD48279 is much larger than its car- 
bon depletion (see Table |4] and above). This is not what we see. 
We thus conclude that the wind of HD48279 is weaker than that 
of HD46966. The value of the mass loss rate is low enough that 
the Ha profiles are photospheric in both stars, and thus simi- 
lar since T^s and \ogg are the same. HD46966 has a luminosity 
L = lO^^^L© while for HD48279, L = lO'^ '^^Lo. It is thus nat- 
ural to have a weaker wind in HD48279. This is what both the 
UV diagnostics and Ha indicate. We note however the difference 
in mass loss rate between the two stars is 0.8 dex if the former 
are used, while it is only 0.4 dex in the case of Ha. This raises 
once again the question of the origin of the differences between 
the UV and optical wind sensitive lines. Whatever the reason, 
we see that two stars with very similar effective temperature and 
gravity can have different wind properties with only a difference 
in luminosity by less than a factor of two. 

5.3. Chemical evolution and rotation 

Rotation is known to affect the angular momentum and chemical 
elements transport through mixi ng processes (meridional circu- 
lation and shear turbulence, see 'Mae der & Mevne3l2000 >). The 
first analysis of surface CNO abundances of O stars and B su- 
pergiants tends to confirm qualitatively the theoretical predic- 
tions: more evolved stars have a hig her N content as a result 
of CNO processing aiid mix ing (Hill ier et"aDl2003t iHeap et aP 
l2006l ICrowflieretal] l200g). However, the role of rotation in 
explaining the chemical patterns of B stars in the Magellanic 
Clouds has been questioned by the results of the first ESO/VUT 
Large Progra mme dedicate d to massive stars (* Evans et"al1l2005L 
2006). Hunte r et al.l (l2008h showed that 60% of their sample of 
B main-sequence stars in the Large Magellanic Cloud displayed 
nitrogen surface abundances expected by evolutionary models 
including rotational mixing, but the remaining 40% escaped 
such an agreement. These outliers are divided in two groups with 
approximately the same number of members: one with highly 
enriched but slowly rotating stars; the other with barely enriched 
fast rotators. The fo rmer group was also identifie d in the SMC 
(Hu nter et alj|2009l) . In the Galaxy, die results of iHunter et al.l 
(200W point to a lack of surface nitrogen enrichment in all main- 
sequence objects, probabl y due to the r ela tively unevolved status 
of most stars. In contrast. iMorel et all (I2006ll2008l) reported the 
detection of N-rich non-supergiant Galactic B stars with low ro- 
tation rates, lending support to the idea that rotational mixing 
cannot explain the surface abundances of a significant number 
of massive st ars also in the Galaxy . 

However, iMaeder et all (l2009l) argued that surface chemical 
enrichment is a multivariate function and does not depend only 
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Fig. 4. Comparison between the best fit models of the UV spectrum (red lines) and models with higher mass loss rates (blue Unes). 
The solid lines are the observed spectra. The left panel shows the comparison for HD 46150, the right panel for HD 46202. The 
narrow absorption features in the C iv 1548-1550 profiles are of interstellar nature. 



on projected rotational velocity. The higher the mass of a star, 
the larger the nitrogen enrichment. Similarly, lower metallicity 
stars have larger surface enrichment at a given evolutionary state. 
In addition, more evolved stars show larger N/H ratios. Hence, 
a combination of ages, masses, metallicities and rotational ve- 
locities can lead to various surface enrichments. It is important 
to disentangl e the different effects to draw reliable conclusions. 
According to lMaeder et alj (f2009). if this is done, rotational mix- 
ing explains most of the observe d abundance patterns. 20% of 
the sample of Hunter et al.' ('2008') might escape this scenario due 
to binary evolution. Brott et al. (201 lb) built population synthe- 
sis models to take all these effects into account. They investi- 
gated the properties of a population of 7-60 Mq stars (with a 
Salpeter IMF) formed continuously for the last 30 Myr and with 
rotational velocities drawn from the distribution of Hunter et al.l 
(2008). They confirm that about 35% of the stars observed by 
iHunter et al.l (2008) cannot be explained by single star evolu- 
tionary models with rotation. But the remaing - and majority - 
of the p opulation is we l l acc ounted for by rotating evolutionary 
models. IPrzvbilla et al.l (1201 0) also confirm the role of rotational 
mixing by showing a tight correlation between N/C and N/O in 
main sequence and supergiants stars. However, the theoretical 
enrichments are too large compared to observations. 



Most of the results summarized above have been obtained 
from the analysis of surface abundances of B main-sequence 
stars. In this paper, we present results for main- sequence O 
st ars. Our sample is smaller than those of iMorel et al.l (l2008h 
or lHunter et al.l (l2008h but it allows to identify some trends, and 
extend their work to the most massive stars. Fig. |6] shows the 
nitrogen surface abundance as a function of rotational veloc- 
ity. The results of Morel et al. (2008), Hunter etal. (2008) and 
iPrzvbiUa et all (l2010h have been added for comparison (small 
symbols). The following comments can be made: 



• up to about 120 km s"' a possible trend of larger N/H with 
larger V sin/ is emerging among O stars, although a signif- 
icant scatter exists and more stars are definitely needed to 
draw firm conclusions. We note that our sample is relatively 
homogeneous in age and metallicity. On the other hand it 
covers a relatively wide mass range (15-60 Mq). Both bina- 
ries and single stars seem to follow this trend. 

• the components of the binaries studied in this paper have 
surface N/H values consistent with those of single stars. 
This conclusion applies only to the binaries studied here and 
should not be generalized. 

• the two fastest rotators show very little enrichment if any. 

• on the main-sequence, in the Galaxy, O stars show on aver- 
age larger values of N/H than B stars. 

These results confirm some of the expectations of stellar evolu- 
tion with rotation. First, more massive stars display larger sur- 
face nitrogen enrichment. This is even more clearly illustrated 
on Fig.|7]where the ratio N/H is shown as a function of luminos- 
ity. A clear trend of higher N content for more luminous stars 
is observed. At a given age, this is what evolutionary models 
predict. The dotted lines show isochrones for 2 and 4 Myr, typi- 
cal of the age of our sample stars. They indicate faster evolution 
and consequently faster chemical enrichment for more luminous 
and thus more massive stars. Interestingly, the two fast rotators 
do not appear as outliers in this plot, but fall perfectly in the se- 
quence followed by most of the stars. This means that given their 
mass, these objects are too young to have experienced a strong 
enrichment in spite of their large rotationa l velocity. This con- 
firms the conclusion of iMaeder et al.l (12009 ) that it is important 
to disentangle the various effects at work in the surface enrich- 
ment of massive stars to pin down the true role of rotation. 

Our sample includes binaries. In Fig. |6] and |7] they do not 
show any peculiar property in terms of chemical enrichment 
compared to single O stars. This does not mean that binarity 
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Fig. 5. Comparison between the observed spectra of HD46966 (blue) and HD48279 (red). The spectra of HD46966 have been 
convolved to the same rotational velocity as HD48279. The similar effective temperature and gravity in both stars is obvious from 
the He i, He ii and Balmer lines. On the contrary, HD46966 presents a stronger C iv 1548-1550 line (but a similar Ho- profile) and a 
lower N/C ratio as clearly seen from the C iii 4070, C iii 4634-4640, C iii 5696, C iv 5801-5812 and N iii 4500-4520 lines. 



has no effect on the history of surface enrichment, but that it 
depends on the binary properties. In particular, not all bina- 
ries have experienced mass transfer which would affect both 
the rotational speed and surface composition dLinder et al.l2008t 
[de Mink et al. 2009). With the exception of HD 48099, all the 
binaries of o ur sample have orbital periods larger than 3 days. 
According to Ide Mink et all it takes more than 8 Myr 

for Roche-lobe overflow to occur in such systems. This result 



is valid for binaries of about 20 Mq and obviously depends on 
masses and mass ratios. The binaries of our sample are thus most 
likely too young and/or too separated to have experienced signif- 
icant mass transfer This conclusion only applies to the systems 
included in our sample and should not be generalized to all bi- 
nary stars. Once mass transfer has occured, it is likely that addi- 
tional mixing processes take place and affect the surface compo- 
sition. Further analysis of such interacting systems should shed 
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Fig. 6. Nitrogen surface abundance (in units of 12+log(N/H)) 
as a function of projected rotational velocity. Single O stars 
are the blue squares and binaries the red triangles. The small 
bl ack triangles (pentagons, open circles) a re the Galactic B stars 
of iHunter et all (l2009h (iMorel et al.ll2008l iPrzvbillaet al"]l20Tol) 
with \ogg larger than 3.5 (i.e. comparable to our sample). The 
grey area coiTesponds to the range of solar values from studies 
of iJie Sun, B stars and HII regions (see Table 3 of iHunter et al.l 
I2OO9I) . 



more light on the exact role of binarity on the chemical appear- 
ance of O stars. 

In Fig.|7] the only clear outlier is star HD48279 with a very 
strong enrichment in spite of a moderate luminosity. To a lesser 
extent, HD46573 can also be considered rather enriched com- 
pared to the expectations of single star evolutionary models due 
to the upper limit on its luminosity. For HD48279 it does not 
seem possible that rotation only can expla i n the large N con- 
tent. According to Fig. 1 of iMaeder et al.l (l2009l) . an increase 
of the rotational velocity from 150 km s"' to 300 km s"' leads 
to an extra enrichment of the order 0.5 dex. HD48279 is 0.7 
dex more N rich than the other stars of the same luminosity i n 
Fig. Q Binarity is excluded for that object ( Mahv et aLll2009h . 
The only possible explanation to the strong N enrichment is that 
HD48279 was member of a binary system in the past and ex- 
perienced mass transfer and/or tidal interaction. Either the bi- 
nary has been disrupted (by supernova kick for instance) or the 
companion is a low-mass compact object rendering any radial 
velocity variation undetectable (see de Mink et al. 2011). In the 
case of HD46573, a large initial velocity or binarity can explain 
the possible slightly larger N/H ratio compared to stars of the 
sar ne luminosity ran ge (HD46573 has been classified as SBl 
bv lMahvetd]l2009b . Finally, HD46150 also only has an up- 
per limit on its luminosity. But even accounting for an extreme 
reduction by 0.3 dex, the N/H content would still be consistent 
with single star evolutionary tracks. 



Fig. 7. Nitrogen surface abundance (in units of 12-i-log(N/H)) as 
a function of luminosity. Single O stars are the blue squares and 
binaries the red triangles. Evolutionary tracks - shown by the 
solid, dot-long dashed, d ot-dashed, dashed lines fo r M=15, 20, 
25 and 40 Mq - are from lMevnet & Maedeii (l2003h . The dotted 
lines are isochrones for 2 and 4 Myr. 
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Fig. 8. Ratio of spectroscopic to evolutionary masses versus 
spectroscopic masses for the sample stars. The dotted line is the 
one to one relation. Black circles are single stars, red squares are 
members of binary systems. The insert shows the position of the 
primary component of HD48099. 
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5.4. Stellar masses 

The determination of stellar masses i s a long-standing problem . 
The so-c alled "mass discrepancy" (iGroenewe gen et aTl Il989t 
iHerrero et al. 1992) refers to the systematic overestimate of evo- 
lutionary masses compared to spectroscopic masses. The for- 
mer are masses derived from the HR diagram and evolutionary 
tracks. The latter are the masses resulting from the knowledge of 
the gravity and radius obtained from spectroscopic analysis with 
atmosphere mod els. The revision o f the effective temperature 
scal e of O stars dMartins et al.ll2002l) has lead to improvements 
(see lGiesll2003l) . Howev er, differences of up to a factor of three 
stm ex ist in some cases. iMassev et aP (l2005h and lMokiem et al.l 
ilOO 'T) reported a mass discrepancy for LMC stars. The latter 
indicate that part of the discrepancy could be caused by fast ro- 
tation in He-rich stars. Fast rotating stars tend to evolve towards 
larger luminosities for the same initial mass as non rotating stars. 
Consequently, using non rotating tracks to analyze them leads to 
overestimates of their evolutionary masses. However, even non 
enriched objects (hence not likely to rotate rapidly) show the 
mass discrepancy. 

In Fig. [8] we show the ratio between spectroscopic (Mj^ec) 
and evolutionary (M^,,,) masses for our sample of Galactic stars. 
Black circles are single stars and red squares the components of 
binary systems. The spectroscopic masses have been computed 
from the effective gravity (that we determine through line fitting) 
corrected from the effect of centrifugal forces. However, because 
of the high uncertainty on the gravity of the SB2 systems, we do 
not include the components of HD 46149 and HD48099 in this 
discussion. In spite of the small number of object, the trend for 
a mass discrepancy at low masses (M < 25 Mq) is present. At 
higher masses, a good agreement is observed, but the eiTor bars 
are large. These resul ts are qualitatively consistent with those of 
iRepolust et al.l (l2004l) who showed that stars with masses lower 
than about 50 Mq seemed to follow a relation parallel to the 1 : 1 
relation in the Mg,, versus Ms„e c diagra m. 

Recently, Weidner & \^ (l2010l) conducted a study of a 
sample of Galactic O stars and concluded that the mass discrep- 
ancy was solved. In their analysis, they mainly focused on bi- 
nary components. They established the evolutionary masses of 
their stars by usin g a spectral type - mass relation. They used 
the Teff-scale of Martins et al. (2005a) to determine an effective 
temperature and luminosi ty for each star and subse quently used 
the evolutionary tracks of lMevnet & Maede^ (l2003h to assign an 
evolutionary mass. They then compared these estimates to dy- 
namical masses resulting from fits to velocity and light curves 
of binary systems. They found a good agreement between both 
types of masses (see their Fig. 3). However, the comparison with 
spectroscopic masses (not shown directly by Weidner & Vink 
but available from their Table 3) still reveals that evolutionary 
masses are larger. Hence, it seems that at least part of the origin 
of the discrepancy is rooted in an underestimate of gravities. 

6. Conclusion 

We have analyzed ten OB stars (nine O stars and one B star) 
in the NGC2244 cluster and MonOB2 association, including 
two bona-fide binaries and one candidate. Optical and UV spec- 
troscopy have been obtained. Atmosphere models computed 
with the code CMFGEN have been used to derive the main stel- 
lar and wind parameters of the target stars. The main results are 
the following: 

• all stars have an age between 1 and 5 Myr, with a trend of 
lower age for the most massive stars. Whether this is a real 



effect or an artifact of the method used to derive the stellar 
parameters is not clear at present. 

• we confirm the existence of weak winds in the latest type 
stars of our sample. For those stars, UV mass loss rates are 
lower than the single line driving limit. UV mass loss rates 
are systematically lower than Ha mass loss rates. Vorosity 
might be a solution to these issues, although its effects on 
dynamics are still unclear. 

• nitrogen surface abundances in O stars indicate a larger range 
of enrichment compared to B stars. A clear trend of larger en- 
richment in higher mass stars is observed, consistent with the 
prediction of evolutionary models for a population of stars 
of the same age. There is no clear difference between single 
and binary stars in our sample, most likely because the latter 
have not (yet) experienced mass transfer or too severe tidal 
effects. Fast rotating non enriched O stars are understood as 
too young to have had time to bring enough fresh nitrogen to 
their surface. 

• the so-called "mass discrepancy" still exists for the latest 
type O dwarfs. Above 25 Mq, evolutionary and spectro- 
scopic masses are in agreement, within the uncertainties. Part 
of the discrepancy seems related to underestimates of the 
spectroscopic masses. 

We need to analyze more stars of the same age to confirm the 
trends observed, especially concerning the surface nitrogen en- 
richment. Winds of massive stars are complex, and a better un- 
derstanding of clumping, both observationally and theoretically, 
is necessary to quantify their properties. 
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Appendix A: Best fit to the observed spectra. 

Figures [ATI to I A.4l show the observed spectra in black together 
with the best fit CMFGEN models in red. 



16 



F. Martins et al.: A quantitative study of O stars in NGC2244 and the Mon OB2 association 



1 

0.8 
0.8 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 




1 . . . . 






1 

0.8 



0.9 
0.8 
0.7 



4800 4850 4900 4950 5000 5050 



1.1 
1 

0.9 
0.8 
0.7 



1 — 1 — ■ — r 



T 



T 




5600 



(a) HD 46573 



1.1 
1 

0.9 
0.8 |- 
0.7 p- 
0.6 



2. 0.0 

1 1.1 
^ 1 

0.9 




0.8 E- 
O.V =- 




1 

0.9 
0.8 
0.7 



: 1 : 


- . 1 





1 


. i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i 




0.95 










0.9 











5600 

x[A] 

(b) HD 46485 



Fig. A.l. Best CMFGEN fits (red) of the optical spectra (black) of stars HD46573 (top), HD46485 (bottom). The models have the 
mass loss rate best reproducing the UV features. The unfitted feature around 4425 A is a DIB. 
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Fig. A.2. Best CMFGEN fits (red) of the optical spectra (black) of stars HD46056 (top), HD46966 (bottom). The models have the 
mass loss rate best reproducing the UV features. The unfitted feature around 4425 A is a DIB. 
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Fig. A.3. Best CMFGEN fits (red) of the optical spectra (black) of stars HD46223 (top), HD46202 (bottom). The models have the 
mass loss rate best reproducing the UV features. The unfitted feature around 4425 A is a DIB. 
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Fig.A.4. Best CMFGEN fits (red) of the optical spectra (black) of the primary (top) and secondary (bottom) components of 
HD46149. The models have the mass loss rate best reproducing the UV features. The unfitted feature around 4425 A is a DIB. 
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